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Observation of B + — > p + p° 
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We report the first observation of the charmless vector- vector decay process B + — > p + p°. The 
measurement uses a 78 fb" 1 data sample collected with the Belle detector at the KEKB asymmetric 
e + e~ collider operating at the T(4S) resonance. We obtain a branching fraction of B(B + — > 
p + p°) — (31.7 ± 7.1(stat.)^g'7(sys.)) x 10 -6 . An analysis of the p helicity- angle distributions gives 
a longitudinal polarization of F L /T = (94.8 ± 10.6(stat.) ± 2.1(sys.))%. 

PACS numbers: 13.25.Hw, 14.40.Nd 



Charmless B meson decays to two pseudoscalar mesons 
(PP) or to pseudoscalar plus vector meson (PV) final 
states have been studied in some detail Q. However, 
measurements of decays to charmless vector- vector (VV) 
final states are rather limited; to date only B — > <j>K* 
decays have been observed 0- The VV decays provide 
opportunities to search for direct-CP and/or T violations 
through angular correlations between the vector meson 
decay final states 03- In B + — > p + p° decays, isospin- 
breaking processes such as electroweak penguins 5] or p°- 
oj interference 0, Q , which may produce a sizable direct- 
CP- violating asymmetry, are expected to be enhanced 
relative to CP-conserving processes such as gluonic pen- 
guins, which are nominally forbidden by isospin symme- 
try. The branching fraction for this process is predicted 
to be 0(1(T 5 ) 0,0. 

In this paper, we present the first observation of the 
VV decay mode B+ -> p+p 0. These decays produce 
final states where both p mesons are either longitudinally 
or transversely polarized. We denote the corresponding 
amplitudes as Hqq and -ffn, respectively. 

The analysis is based on a 78 fb" 1 data sample contain- 
ing 85 x 10 6 B meson pairs collected at the T(4S*) reso- 
nance with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider. We also use an off-resonance data 
sample of 8.3 fb _1 collected at a center-of-mass energy 
that is 60 MeV below the T(4S) resonance. 

The Belle detector is a large-solid-angle magnetic spec- 



trometer. Charged particle tracking is provided by a 
three-layer silicon vertex detector and a 50-layer central 
drift chamber (CDC). Charged hadron identification is 
provided by dE/dx measurements in the CDC, and ar- 
rays of aerogel threshold Cerenkov counters (ACC) and 
time-of-fiight scintillation counters (TOF) that surround 
the CDC. An electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) provides photon detection and 
electron identification. All of these devices are located 
inside a superconducting solenoidal coil that provides a 
1.5 T magnetic field. An iron flux-return located out- 
side of the coil is instrumented to detect K® mesons and 
muons. The detector is described in detail elsewhere 

We select B + — > p + p° candidate events by combining 
three charged pions and one neutral pion. We require 
that each charged track has a transverse momentum p t > 
0.1 GeV/c, and is consistent with originating from within 
Sr < 0.1 cm in the radial direction and \Sz\ < 5 cm in 
the electron beam direction of the run-by-run-determined 
interaction point. We also require that the three charged 
tracks be positively identified as pions by the CDC, ACC 
and TOF systems. 

Candidate tt° mesons are reconstructed from pairs 
of photons with an invariant mass in the range 
0.118 GcV/c 2 < M(77) < 0.150 GeV/c 2 . For the ECL 
barrel region (32.2° < 9 < 128.7°), photon energies 
greater than 50 MeV are required; for the ECL endcap 
region (17.0° < 6 < 31.4° or 130.7° < 6 < 150.0°), this 
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requirement is increased to 100 MeV. In addition, we only 
accept 7T° candidates with a T(45) center-of-mass system 
(cms) momentum > 0.5GeV/c. The 7r° candidates 
are kinematically constrained to the nominal 7r° mass. 

Candidate p mesons are reconstructed via their p° — > 
7r + 7r _ and p + — » tt + tt° decays. For both the charged 
and neutral modes, we require 0.65GeV/c 2 < M{-ktt) < 
0.89GeV/c 2 . 

B + — > p + p a decays arc identified using the beam- 
energy constrained mass M& c = a/ (-Ebcam) 2 — (pb) 2 and 
the energy difference AE = Eb — -Ebcam, where -Eboam 
is the cms beam energy, and ps and Eb are the cms 
momentum and energy, respectively, of the B + — > p + p° 
candidates. The AE distribution has a tail on the 
lower side caused by incomplete longitudinal contain- 
ment of electromagnetic showers in the Csl crystals, and 
the AE resolution varies slightly depending on the 7r° 
momentum. We select events in the region \AE\ < 
0.4 GeV, Mbc > 5.2 GeV/c 2 with a signal region defined 
as -0.10 GeV < AE < 0.06 GeV and 5.272 GeV/c 2 < 
Mbc < 5.290 GeV/c 2 . The p — > -kit daughter momen- 
tum distributions are helicity-state dependent: for the 
Hqq state, the two pions have an asymmetric momentum 
distribution where one pion has low momentum (in the 
range ~ 1.3GeV/c) while the other has high momen- 
tum (1.3 ~ 2.8GeV/c); for the Hi\ state, the two pions 
tend to have the same momentum. Because of its higher 
probability for having a low energy pion, the Hqq state 
has a lower reconstruction efficiency and a AE resolution 
that is, on average, about 15% broader than that for the 
H\x state. 

There are large backgrounds from e + e~ — ► qq con- 
tinuum events (q — u,d,s,c), which tend to have a 
two-jet-like structure. These are suppressed by requir- 
ing | cos#thr| < 0.8, where t hr is the angle between the 
thrust axis of the candidate tracks and that of the re- 
maining tracks in the event. We achieve further suppres- 
sion by means of a likelihood ratio derived from a Fisher 
discriminant formed from six modified Fox- Wolfram mo- 
ments [ll| and 9b, the angle between the B flight di- 
rection and the electron beam direction. The combined 
rejection for continuum events is 98%, with a 65% loss in 
signal. 

Background contributions from b — > c processes are 
investigated with a large sample of Monte Carlo (MC) 
events, for which no signal-like peak is found in either 
the AE or Mf, c distributions. Some rare B decay pro- 
cesses, such as B + — > i]'p + , K* + p°, p + K*° and pir, can 
survive the event selection but are displaced from the 
signal in AE. Moreover, these modes have small branch- 
ing fractions 0] and low reconstruction efficiencies. MC 
estimates based on measured upper limits for the branch- 
ing fractions indicate a possible signal-region yield from 
these rare modes of seven events; this is taken into ac- 
count in the systematic error determination, as discussed 
below. 



Figure^(left) shows the AE projection of the selected 
events in the 5.272 GeV/c 2 < M bc < 5.290 GeV/c 2 sig- 
nal region. The curve shows the results of a binned 
maximum-likelihood fit with three components: signal, 
continuum background, and BB background. The signal 
is represented by the sum of a Gaussian and a "Crystal 
Ball" line shape (CB) function with parameters de- 
termined from an Hqq signal MC that is calibrated with 
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"7r° events. A linear function 



with a slope determined from the off-resonance data is 
used to represent the continuum background. The BB 
background contribution is modeled by a smoothed his- 
togram with a shape that is obtained from MC. In the fit, 
all parameters other than the normalizations are fixed. 

The fit gives a signal yield of 58.7 ± 13.2 events. 
The statistical significance of the signal, defined as 
y/—2hi(Co/C max ), where £ max is the likelihood value at 
the best-fit signal yield and Cq is the value with the signal 
yield fixed to zero, is 5.3. 
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FIG. 1: AE (left) and M hc (right) fits to the B + -> p+p° 
candidate events. The signal component is shown as a dot- 
dashed line. The sum of BB and continuum components is 
shown as a dashed line. The shaded histograms represent the 
BB background. 

Figure n (right) shows the Mf, c projection of events in 
the -0.10 GeV < AE < 0.06 GeV signal region. The 
curve shows the results of a binned maximum-likelihood 
fit that uses a single Gaussian with a MC-determined 
width to represent the signal, a threshold (ARGUS) func- 
tion [lij l for the continuum background with shape pa- 
rameters that are determined from the AE sideband (de- 
fined as 0.1 GeV < AE < 0.4 GeV), and a smoothed 
histogram obtained from MC to represent the BB back- 
ground, normalized according to the MC expectation. 
This fit gives a signal yield of 49.1 ± 9.7, with a statistical 
significance of 6.5. 

The fit results and MC efficiencies are summarized in 
Table [I] Since the AE distribution provides stronger dis- 
crimination against rare B-meson decay backgrounds, we 
use the AE fit result for the branching fraction calcula- 
tion. 

Figure [21 shows signal yields extracted from fits to the 
AE distributions for different A/(7r+7r~) and M(7r+7r°) 
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TABLE I: Signal yields from the fits to the AE and M bc 
distributions together with the MC-determined efficiencies. 





AE fit 


M bc fit 


Yield 


58.7 ± 13.2 (5.3ct) 


49.1 ±9.7 (6.5a) 


Efficiency eoo 


2.11% 


1.59% 


Efficiency en 


3.45% 


3.07% 



mass bins; the irir mass spectra from the signal MC are 
shown as shaded histograms. The data agree reasonably 
well with B + — > p + p° MC expectations. 
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FIG. 2: Data points are the results of fits to the AE distribu- 
tions for each M(ji + tt~) bin (left) and M(7r + 7r°) bin (right), 
where D — > K-k, tvtt, tttv events are excluded by —120 < 
M{wk) - M D o < 30 MeV/c 2 and |A/(7T7r°) - M D +\ < 50 
MeV/c 2 vetoes; the histograms are expectations from the sig- 
nal MC. 

We examined the possible contribution from 
non-resonant processes using MC-generated 
B + — > 7r + 7T _ 7r + 7r° events where the final states are 
distributed uniformly over phase space. After the 
application of all selection requirements, including the 
p mass cuts, we find an efficiency that is less than 2% 
of that for B + — » p + p° decays. Possible contributions 
from B + — > a^7r° or B + — > aj ) 7r + decay are examined 
and found to be smaller than those from non-resonant 
decays. To account for these contributions, we perform 
X 2 fits to the distributions shown in Fig. [2] with a p plus 
non-resonant 7T7t component included. The resulting 
non-resonant yield increased by la is included in the 
systematic error. 

We use the p — > 7T7T helicity-angle (6*hei) distributions 
to determine the relative strengths of Hqo and Bl\\. Here 
#hci is the angle between the p flight direction in the B 
rest frame and tt + flight direction in the p rest frame. The 
signal yields determined from fits to the AE distributions 
for each helicity-angle bin are plotted versus cos#hei in 
Fig.|3]for the p° (left) and the p + (right). In determining 
these yields, bin-by-bin 7?oo-signal-MC-determined shape 
parameters are used to extract the signal. We perform 
a simultaneous x 2 fit to the two background-subtracted 
p helicity-angle distributions using MC-determined ex- 
pectations for the Hqq and Hn helicity states. The 



fit results, shown as histograms in Fig. [3J indicate that 
the longitudinal (Hqo) state dominates. We obtain the 
acceptance-corrected longitudinal polarization ratio 

-± = (94.8 ± 10.6(stat.) ± 2.1(sys.))%, 

where the systematic error includes uncertainties in the 
signal yield extraction and the polarization-dependence 
of the detection efficiency This dominance of -f/oo is con- 
sistent with theoretical predictions [lfij . 
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FIG. 3: Data points show the background-subtracted cos #hci 
distributions for the p° (left) and p + (right). In each plot the 
dashed (dot-dashed) histogram is the Hoo (Hn) component 
of the fit; the solid histogram is their sum. The low yield of 
events near cos# hol ( p +) = 1 is due to the p x o > 0.5GeV/c 
requirement. 

We use MC-determined efficiencies for the two helicity 
states and the measured polarization ratio to determine 
the branching fraction. We assign a 3.4% systematic er- 
ror for tracking that is obtained from a study of partially 
reconstructed D* decays; a 3.6% error for the particle 
identification efficiency that is based on a study of kinc- 
matically selected D* + — > D°ir + , D° — > K~ix + decays; 
a 4.0% 7T° reconstruction systematic error that is deter- 
mined from comparisons of rj — > 7r°7r 7r° to rj — > 7r + 7r~7r° 
and rj — ► 77; a 5.4% error for continuum suppres- 
sion that is estimated from a study of B + — > D°ir + , 
D° —> K + ir~Tr° decays; an error associated with the AE 
fit of tg'7% that is obtained from changes that occur 
when each parameter of the fitting functions is varied by 
±1<t; a 1% error for the uncertainty in the number of BB 
events in the data sample; a i? 6 7 % error to account for 
a possible contribution from non-resonant decays; and 
a 3.3% error due to uncertainties in the rare B decay 
background that is estimated from the change produced 
by fitting the AE distribution with the inclusion of an 
additional component normalized at its MC expectation. 
We also include a error due to the uncertainty in 

the mixture of helicity states that is obtained from the 
changes in the branching fraction that occur when the 
longitudinal polarization ratio is lowered by la or raised 
to 100%. The quadratic sum of all of these errors is taken 
as the total systematic error. We obtain the branching 
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fraction 

B{B+ -^p + p°) = (31.7±7.1(stat.)l|;?(sys.)) x 10~ 6 . 

As a check, we examined the decay mode B + — > p + D°, 
D° — » ir + ir~, which has the same final state particles 
as the mode under study, including a 7r° with a similar 
momentum distribution. The same analysis procedure 
is applied except for an \M(irir) — M D o\ < 13 MeV/c 2 
mass selection. For this mode, we obtain a signal yield of 
42.3 ±8.5 events, consistent with MC expectations based 
on the known branching fraction values Q . 

Direct CP violation would be indicated by a difference 
in partial rates for B~ — > p~ p° and B + — > p + p°. Sep- 
arate fits to the AE distributions find 29.3 ± 9.5 p~ p° 
events and 29.3 ±9.1 p + p° events. Since backgrounds 
from generic BB decays should contribute equally to 
p~ p° and p + p°, we fix the normalizations for BB com- 
ponents at half the value determined from the combined 
fit. 

The charge symmetry of the detector and reconstruc- 
tion procedure is verified with a sample of B + — > D°ir + , 
D° — > K + -k~-k° decays and their charge conjugates. 
Here the analysis procedure is similar to that for B + — » 
p + p° but replacing one ir + by a K + and the invariant 
mass requirements by \M(Kmr°) — M D o\ < 50MeV/c 2 . 
For these events we find a partial rate asymmetry of 
(—2.1 ± 2.5)%, which is consistent with zero. We as- 
sign 2.5% as the systematic error for the detection and 
reconstruction asymmetry. The systematic error associ- 
ated with the AE fitting procedure is determined to be 
(-1 2)% by shifting each parameter of the fitting func- 
tions by ±ler and taking the quadratic sum of the result- 
ing changes in Acp ■ The quadratic sum of these errors 
is taken as the total systematic error. We obtain the CP 
asymmetry 

A (m , „ T „(n - N (P-P°) ~ N (P+P°) 

A C p{B+ -> p+p ) = — — — 

= 0.00 ± 0.22(stat.)±0.03(sys.). 

In summary, we have observed the decay B + — > p + p° 
with a statistical significance of 5.3. We measure 
the branching fraction to be B(B + — * p + p a ) = 
(31.7±7.1(stat.)lg 7(sys.)) x 10" 6 , where the systematic 
error includes the error associated with the helicity- 
mix uncertainty. An analysis of the helicity-angle 
distributions gives the longitudinal polarization ratio 
T L /T = (94.8 ± 10.6(stat.) ± 2.1(sys.))%. We also mea- 
sure A C p(BT -> pT p 0) = 0.00±0.22(stat.)±0.03(sys.). 
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